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Abstract 

The relativistic model with two types of planar fermions interacting with the 
Chern-Simons and Maxwell fields is applied to the study of anyon superconductor. It 
is demonstrated, that the Meissner effect can be realized in the case of the simultane- 
ous presence of the fermions with a different magnetic moment interactions. Under 
the certain conditions there occures an extra plateau at the magnetization curve. In 
the order under consideration the spectrum of the electromagnetic field excitations 
contains the long-range interaction and one massive "photon" state. 



1. Introduction 

The zero-temperature Meissner effect presented in the 2+1 dimensional anyon 
matter, provoked a considerable efforts in order to promote the Chern-Simons gauge 
theory as a hypothetical candidate for the high T c superconductivity. 

The most important points in that development are existence of the massless 
(Goldstone) pole in the current correlators |], cancelation of bare and induced C-S 
terms ||, and detailed calculations of effective action and thermodynamical potential 
for the fermions ineracting with C-S and Maxwell fields J| f|, [|, ||. 

Among the others it was shown, that the Meissner effect is partial, i.e. magnetic 
field starts to penetrate into the sample at any non-zero temperature ||. The second 
result concernes the effective mixing of Maxwell and C-S fields and occurence of two 
distinct massive excitations [f|. 

In the present paper we try to give some complementary insights into those 
intriguing questions. 

It is demonstrated, that the Meissner effect exists only if the matter consists 
of the two types of fermions with opposite signs of magnetic moment interaction. 
Such a system can be naturally realized considering the planar relativistic fermions. 
Note, that the different versions of the relativistic anyon superconductivity have been 
considered in || [/]. || . 

Studying the gauge field propagators, we found out, that results cited in can 
be understood as a consequenses of momentum expansion of the certain structure 
functions. Rising the order of terms taken into account, it can be shown, that there 
exists the single massive excitation. 

The paper is organized in the following way. Section 2 containes the description 
of model and some basic definitions. 

In Section 3 we present the thermodynamical potential and relevant physical 
quantities. 

In Section 4 we show, that the Meissner effect can take place if the system contains 
two types of fermions. Further we calculate the dependence of critical magnetic 
field on the temperature and establish the existence of additional plateaux at the 
magnetization curves. 

In Section 5 we derive the electromagnetic field propagators and analyze the 
spectrum. 



2. The Model 

Consider the system of Dirac fermions in 2+1 dimensions interacting with Maxwell 
electromagnetic field and Chern-Simons gauge field. 



1 



The total lagrangian for this system is a sum of matter and gauge field lagrangians 
4/i ^ h 2tt 



C = - — F IW F I1V ---^ e^\d u a x + ecn e A + 



+ ^{ihcj^D^ — amc 2 }^, 

D li = d ft + i(e/h)(A fl + a, i ), m > 0, a = ±1, 

where 5 ~ lCT 9 m is the interplanar distance in the multilayered system and its 
presence is justified by the use of the three dimensional Maxwell field and electric 
charge. The Chern-Simons field is in fact two-dimensional, but will be measured in 
three dimensional units, ip is the 2-component fermion field and en e is the planar 
density of the background neutralizing charges. In 2+1 dimensions the sign of the 
mass term is essential and describes the helicity of the fermion. 

In what follows, we study the static properties of this relativistic system. Analysis 
will be performed in the self consistent field approximation, developed in ||. It means, 
that the gauge fields must satisfy the equations 

- —AA + ecn e = (J°(r)), (1) 

- -e kn d n B = (J k (r)), (2) 
A*o 

£^6=<j0(r)>, (3) 
n 7r 

f ^^ n a =(J fe (r)>, (4) 
where the thermal averages are defined in terms of the grand canonical ensamble 

(• • •) Tr e -(^e-Atmc2JV)/A B T _ Tr | _ _ e -(H e -^mc' 2 N)/k B T^ 

Here ks, T and fi are the Boltzmann constant, temperature and the dimensionless 
chemical potential respectively, while A" denotes the particle number operator 

N = J ^ f (r)^(r)dr. 

The fermion hamiltonian is given by 

H e = J ip^ (r){ihcy°jkDk(r) + ecA (r) + eca (r) + amc 2/ y°}ijj(r)dr. 
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Introduce the thermodynamic potential 

^ ,m a \ , ™, m f #e(A a) - umc 2 N) 
Q e (T, /i, A, a) = -k B T In Tr exp I - el 1 * \ , 

in terms of which the current averages are given by 

The gauge fields are presented as a sums of background and fluctuating parts 

A u = A h u + A^, a n = a° u + a^. 

Backgrounds correspond to the uniform Maxwell and Chern-Simons magnetic fields, 
i.e. Aq = Oq = 0, B h = const, b h = const. 

Separate the fermionic hamiltonian into the dimensionless free and interacting 
parts 

H e = mc 2 (H + H int ), 
H = j ^(r){t£ ol ° lk D h k (r) + <T 7 °}V(r)dr, (6) 
D b k (r) = d k + i(e/h)(A h k + a b k ), 



^ = -^/^(r){4(r) + <(r)}dr, 



where £q = h/mc is the Compton wave length for the fermions. 

Applying the perturbation theory formalism, we get for the thermodynamic po- 
tential 



Q e = fi - k B T\n(rexp j-J^ H int {j)dT 



mc 2 



o 



Qo is the thermodynamic potential for the system in the uniform magnetic background 
Q (T, fi, A h , a h ) = -k B T\nTiexp{-(3(H (A h , a h ) - fiN)} (8) 



and (• • -) is defined as 



(• • -) Tr { e -^o-^)J = Tr {• • • e ^ H ^' 
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In (7) H int {r) is the interaction hamiltonian in the Matsubara representation 

and T denotes r-ordering. Matsubara fields are given by 

iKr, r) = eF £{ V „ p (r) e -<^ + ^ np (r)e^}, (9) 

$(t, r) = e"^ E{<^np(r)e^ + b np v np {v)e~^}. (10) 

Remark, that except the lowest energy eigenvalues the solutions of Dirac equation 
are always paired. For the lowest energy there is asymmetry, i.e. there is no vq mode 
for as = +1 and no u mode for as = —1, where e = sgn(ei? b + eb h ). 



3. Thermodynamic potential 

In this item we shall find the analytical expression for the thermodynamic poten- 
tial in the second order approximation with respect to the gauge field fluctuations. 
All the operator expressions such as hamiltonian, currents, etc. are assumed to be 
normal ordered. 

Substituting H into (8) and taking into account the degeneracy of the Landau 
levels we get 

Q (A h ,a h ) k B T ^ + k B T ^ 

— 4^T- = 2^ ? ln(1 " M + 2^ ? ln(1 " P J - (U) 

where p+ and p~ are the Fermi distribution functions for the particles and antiparticles 
respectively 

Pn = {l + exp[/3(w„=F/i)]} _1 , u n = VT+ 2hn, 



and £ is the magnetic length. 

Thermodynamic potential in the second order approximation is given by the 
functional 

n e (A, a) = n (A h , a h ) + ec J IP(r){4(r) + <(r)}dr+ 



+ ^-J n^(r 1? r 2 ){4( ri ) + a^( ri )}{4(r 2 ) + a[ (r 2 )}d ri dr 2 , 
where IT M (r) and IP"(ri, r 2 ) read as 

nM ( r ) = ^2 S^WnTX - 2^ E*VVA» = co ^; 

n^( ri , r 2 ) = iy ^e^-^n^k). 

Polarization operator possesses the following tensor form 

n 00 (k)=^n E (k 2 /2), 

n ofc (k) = ■^ie fc, ife , ncs(k 2 /2), 

U kl {k) = 27ce km e ln k m k n U M {k 2 /2), 

where IIe, lies and I1m are the structure functions, presented in the Appendix. 
Explicit calculations show that 

n »= ^rt-^pEr-n. n t = o (13) 

and the final expression for the thermodynamic potential is given by 
Q e (A, a) = n (A b , a h ) + ecll J {A f (r) + 4(r)}dr+ 



+ 



e 2 m 2 c 2 



2m h 
e 2 2mc 



J {Al(v) + 4(r)}n E {4(r) + 4(r)}dr+ 



p Zinc f *■ 
+ 2m~~h T ~ J {Ai ° {r) + a o(r)}n cs {5 f (r) + b { (r)}dr+ 

+ ^ / {Bf(r) + fef ( r )}n M {5 f (r) + 6 f (r)}dr, 

where fl stands for the differential operator n(— £ 2 A/2) (A is the Laplace operator). 
Using (5), the current averages can be expressed as a linear functions of gauge field 
fluctuations 



(j°(r)) = ecn + ^n E (4 + 4)+ ^flcs(5 f + & f )> (14) 



ft 2 "\ » h 



( J*(r)> = -e kn d n |^ n cs + oj ) + ^ n M (5 f + 6 f )| . (15) 



4. Uniform magnetic field and the Meissner effect 

As a starting point, consider the system in the zeroth order approximation, i.e. 
= = 0. Equations of motion (1) and (3) are reduced to 

n e = Il , (16) 
b h = — n e , (17) 

where il is given by (13). 

Note, that n e is the free fermion density in the sample and (16) serves to define 
the chemical potential \x = /x(T, B h , b h , n e ). 

Equation (17) indicates that the Chern-Simons magnetic background is created 
by the free fermion density. For the typical value n e = 10 18 m -2 we get 

2 . i n 7 

|6 b | = —— Gauss. (18) 
I u o\ 

We shall study the compound system, consisting of two sorts of fermions with 
equal gauge couplings e\ = €2 = e and the different helicities o\ = — a"2, corresponding 
to the different signs of the magnetic moment interaction, which on its turn can be 
associated with spin up and spin down fermions. In that case the r.h.s. of the equation 
(16) gets the contributions, defined by (13) from both sorts of particles. Remark that, 
since the different sorts have the equal charges, corresponding magnetic lengths are 
also equal. Moreover, as it follows from (18), realistic values of B h are small compared 
with those of b h and can be neglected: 

e = sgn(eB h + eb h ) = sgn(e6 b ). 

In other words, e can be considered as independent of B h , and without loss of gener- 
ality we can set <7i£i = — a 2 S2 = 1. Taking into account the contributions from both 
types, equation of motion (16) can be rewritten in the following form 

h 

where v\ and u 2 are the corresponding filling fractions 



U e = 7^2 (Vl + V 2), (19) 



^i = E^i)-E^0. (2°) 



n=0 n=l 
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ra=l n=0 

and fj,i t 2 are the chemical potentials for the types 1 and 2 respectively. In (20) and (21) 
we took into the consideration the spectral asymmetry of the one-particle hamiltonian, 
which is reflected in the absence of n = modes in certain terms. 
Introduce the partial contributions to the particle density 



h 

to 

By means of (19) we can express h in terms of the average filling fraction 



n^=™. (22) 



h=^, u= ^±^. (23) 

v 2 v ' 

Using (22), (23), we express the partial filling fractions in terms of the average one 

via = (i) ~~ (2) " (24) 

Tie + Tie 

Substituting h, from (23) and (22) into (20) and (21) we see that the chemical 
potentials depend on temperature and the average filling fraction v. 

The same time, expression for the magnetic length (12) together with equations 
(17) and (23) yields 

1 e eeB h 

- = - + r, 25 

v u 7rn e h 

which reflects the one to one correspondence between v and B h . Consequently, any 
quantity depending on B h can be also viewed as a function of v, and vice versa. 

The value of the background magnetic field is determined by the external mag- 
netic field B cxt and the magnetization M 



B h = B cxt + M(B h ), M(B h ) = - ^ dJ ^\ (26) 



where T = T\ + JF 2 is the Helmholtz free energy density of the composite system. 
The separate contributions are 

Fi,2(B h ) = +/ii, 2 mc 2 (iV li2 )}, 

where (Ni^) are the thermal averages of the fermion numbers determined by 

- — ^2 E PtM - i E PnM = U e\ 



Area 2*$^"™ 2*%^' 
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(N 2 ) 



h 



h 




Area 



2rr£ 2 n 



n=l 



n=0 



Here the quantities h, /xi and /X2 are the functions of v. In order to exhibit the 
global behaviour of T[y) we have used the numerical methods. Consider the case 
of the equal concentrations rvp = = n e /2. Corresponding results are plotted 
in figure 1. As one can see the free energy density of the composite system has a 
local minima at the integer values of the average filling fraction, however the separate 
contributions of each type of particles exhibit no minima. 

The cusp like structure of the free energy is the manifestation of the Meissner 
effect, when the system tryies to expell the magnetic field from inside the sample ||. 

As we see, in contrast with the previous calculations, the Meissner effect does not 
exist in the single fermion system, but only in the composite one, where the diversity 
in the magnetic moment interaction plays the decisive role. 

Magnetization is expressed with the help of the free energy. The later containes 
the additive contributions from the one-particle state energies. One particle hamilto- 
nian considered in || , takes into account the magnetic moment interaction only with 
Maxwell magnetic field. Further simplification is achieved by taking 



assuming the vanishing magnetic moment interaction. 

In our approximation the magnetic interaction term containes the contribution 
from both Maxwell and Chern-Simons magnetic fields. The nonrelativistic limit of 
one-particle energy spectrum is given by 



Remark, that (27) is in fact the half sum of E n (+) and E n (—). Separate use of 
E n (+) or E n {—) does not lead to the free energy with localized minima and only 
their simultaneous contribution has a cusp like form. 

Equations for the chemical potentials cannot be solved in the global form. How- 
ever, one can find the analytic form of /xi^fV) nearby the integer values of u, where 
the free energy achieves the local minima. This enables to analyze the system in more 
details near these minima, where the Meissner effect just takes the place. Below we 
present these calculations for the asymmetric concentrations i.e. when n]p ^ 

Here we shall deal with the minimum corresponding to v = 1 and for concretness 
set e > and e = 1. In that case we have u = 1 and 




(27) 




h = 7cn e £l(l + a), 



eB 



(28) 



a = 



7rn e h 
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Note that for the characteristic values of the internal magnetic field (B < 200 Gauss) 
we have a < 10~ 5 , (1 + a) -1 = 1 — a and present (24) as 

v\ — (1 + «)(1 — a) = 1 — ai, 
v 2 — (1 - a)(l - a) = 1 - a 2 , 



a = c 



Here a measures the asymmerty in the concentrations of the different types of fermions. 

Further simplifications are due to the fact that in the considered range of tem- 
peratures we have (5 > 3 • 10 7 , and consequently for /xi ;2 > 



^^ = l + e ^+M <e ' 107 ' 



meaning that the main contributions to (20) and (21) come from p+(/Xi,2), forcing z/ 12 
to be positive. With this assumption we take \i\ > and present it as 

A*i = 1 + — g — ^ 

where w\ is to be found. Due to (28) one has h ~ 5 • 10~ 7 . Consequently, for the 
Landau levels with 2hn « 1 we can use c<;„ = 1 + hn and write down 



(3{uj n - /^) = (3h (n - ^ + ^) • 



Characteristic values of /3 and /i are such that f3h > 15, allowing to neglect the 
contributions coming from the higher Landau levels. This permits to write (20) as 
u i = Poi.Pi) + ptini) or m the equivalent form 

1 — Oil = 



1 + e-^e^ 1 1 + e^e^ 1 ' 
where 77 = nn e £lP(l + a)/2. From this equation we easily solve w\ as 

e ±nwi = — _ — (Jl + ajsh 2 ^ ± a^hr]) , on = a - a + aa. (29) 

1 =F ai V ' / 



=F on 

The same time, we take 



1 , 3 W2 , 

jl2 — H 
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and performing the same manipulations obtain 

,±rnv2 = — _ — fj x + a 2 sh 2 v ± ) a 2 = a + a - aa. (30) 

F«2 V / 



e 

1=F«2 



So, the leading contributions to the different physical quantities come from the 
following Fermi distribution functions 

1 1 

Po = ; > Pi(Pi) = --> ) (31) 

Pl{P2) = , pt(p2) = , (32) 

where e^ 1 and e^ 2 are defined by (29) and (30). 

In this approximation the system magnetization is M = Mi + M 2 where 



hn + , 



4m<5 {/(« - «) + 1} , (33) 



— PniP2) 



4m<5 {/(« + «)-!}, (34) 



/(*) = 7 In , = = * * e = ™ e ^/3 (35) 

These quantities define the magnetic and thermal properties of the system. 
In the case of the symmetric concentrations (a = 0) one gets 

b [i en e h 1 n _ y/a 2 + 4e~£ — a eB b 



M(T, B b ) = ^ In \ ' -, a 



Behaviour of M as a function of B h is depicted in figure 2 for the different values of 
temperature. Presented curves are the same as given in earlier work ||. 

On figure 3 we present the dependence of B h on B ext which is set by (26). As one 
can see, when the external magnetic field is applyied, the internal one is in general 
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non-zero. In other words, the Meissner effect is not complete. This observation was 
first made in || . The same time the magnitude of B h is quite small until _B ext reaches 
some critical value B cr (T), which depends on the temperature. Above this value 
external magnetic field begins the notable penetration in the sample. 

The critical value of the external magnetic field is evidently related to the small 
interval at B h axis, where the magnetization curve drastically changes its direction, 
i.e. where the curve passes the point of maximal curvature (PMC). The lower is 
temperature, more narrow is the interval and it is easier to establish the corresponding 
critical magnetic field. In order to find the approximated value of B cr (T) we can use 
the behaviour of the derivative 

dM_ _ _ ^e 2 k B T , , h 2 /mkBT yi/2 
dB* ~ rfra ^ > ' 

Until the curve M(B h ) reaches the PMC, its slope can be considered to be con- 
stant and therefore can be determined by its value at the origin, where it is of the 
order of 10 5 for T = 50°K and 10 50 for T = 10°^. On the other hand, in the low 
temperature regime the curve becomes practically horizontal after passing the PMC. 
Obviously, somewhere in the vicinity of PMC one has dM/dB h = — 1 and using this 
relation as the definition of the location of PMC, one gets the corresponding value of 
the internal magnetic field to be 



,1/2 

Hqk 1 ) 




irn e h /rnksT 



Substituting B (T) into B cr = B — M(B ) and keeping the leading terms we get 

2m5 irhd 2TT 2 h 2 n e 5 

Corresponding curve is presented at figure 4. We observe that the value of the critical 
magnetic field decreases for higher temperatures and completely vanishes near T = 
90°K, which seems to be the critical temperature for the system. However, the area 
of the maximal curvature becomes smeared for T = 90°K, and the critical magnetic 
field is ill defined, causing the critical temperature to be also ill defined. On the 
contrary, the PMC is well located for T = 0°K and the corresponding value of the 
critical magnetic field is given by 

B-(0) = 

2m5 

Situation with the asymmetric concentrations is presented at figures 5 and 6. In 
the interval, which includes B h = 0, magnetization practilally vanishes. The width of 
the plateaux, is equal to 2a-(im e h/e) and does not depend on temperature. This kind 
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of behaviour of the magnetization can be explained by the fact, that the magnetization 
curves of each kind of fermions have the standard form set by (35), but shifted by a 
relatively to each other according to (33) and (34). 

The genuine value of the asymmetry parameter a must be determined by the 
equilibrium conditions imposed on the chemical potentials. 



5. Gauge field propagators and mass spectrum 

Equations of motion for the gauge fluctuations can be introduced as the extremals 
of the effective lagrangian 



-±B<B<-^ { B< + b <)n%(B< + b <) + 



+ e ^f^- ^(4 + 4)n-(5 f + 6 f ), 

where the script "tot" means that the corresponding quantities include the contribu- 
tions from both types of fermions. 

In order to study the gauge field propagators we introduce the gauge fixing term 

and present the total lagrangian as 

C cS + C a = - A T DA, A T = (ecA l , eA l n , eca l , ea* n ) . 

Gauge field propagators are defined by the inverse of the matrix D 

DG = GD = 1 
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and can be presented in the following form 



G 



G\ 



oo 



G'l 



reO 



Gf 



Gf 



G\ 



0m 



G™ 



QOm 



Qnm 



Gf 



Gf 



Gf 



Gf 



QOm 



Qnm 



QOm 



Qnm 



Here the tensor blocks are given by 



• c nm l,m 

q On _ _ q nO _ lt ^ q CS 

k 2 



^ni^mj h.i uj h.n urn 

Qnm _ fc fc K ^ (7 M _|_ Q\\ 



k 4 



k 4 



fi e 



7T 



7T 

_ TTtotTTtot 

^0 



^-YCS _ ^ f TTtQt TTtQt TTt0tTTt0t ^0 TTtotl 



i v E 



5k 2 



V/j, e 2 c 2 



IT \ 2 IT 2 

TTtot \ " TTtotTTtot 

77 ii cs J " 772 11 e 11 m 

^0 / ^0 



777 A 

i? = n™n™ - n^re 5 + 



/i e 



+ 



<5 2 k 2 f/x e 2 



7T \ 2 7T 2 

TTtot I / 1 ___ TTtot \ _ TTtotTTtot 

ii M ' I 1 ~~ 77 ii CS J " 772 ii E ii M 
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and (uo/tt) 2 T> is the determinant of D. In the Coulomb gauge (a = 0) we obtain 



G oo (fc) _ 2o E + 3 2 E k 2 + ^k 4 - 



d + rf 2 k 2 + d 4 k 4 + 



Gl n {k) = ie nm k m { ^ + 9 ™ + ^ + ^ + 



k 2 d + d 2 k 2 + d 4 k 4 + 



/ k n k m \ ( n M n M 4- n M k 2 4- n M k 4 4- 



k 2 ; V k rf o + 4k 2 + rf 4 k 4 + ■ ■ ■ J' 

The constants g's are the coefficients of the momentum expansion of the structure 
functions in the nominators of (36), (37) and (38), while ef s are set by D = d 2n k 2n . 
The ratios 

#o + # 2 k 2 + # 4 k 4 + • • • 



d + d 2 k 2 + d 4 k 4 H 



represent the short range interaction. 

The leading contributions to the first two coefficients of the determinant expan- 
sion are given by 



inn 

do = n&(o)ng(o)- Hn^(o), 



mS I 2 dU^\x) 
Hoe 2 2h 2 dx 



x=0 



Using (A.l) - (A. 5), one can show that do and d 2 are the nonvanishing positive 
numbers. Consequently keeping the lowest order terms one can write 

go + g 2 k 2 + g 4 k 4 h g Q g 1 



d + d 2 k 2 + d 4 k 4 + ■ ■ ■ d + d 2 k 2 d 2 k 2 + M 2 ' 

where Ai 2 = do/d 2 > gives rise to the short range interaction. Using the zero 
temperature values of the structure functions (see Appendix), one can estimate the 
corresponding value of the mass parameter 



M 2 = -ZJ2!Z_± = (39) 



h 2 [Xoe 2 n e h 2 
= A 2 ' 



where Al is the London penetration depth. Note, that M. does not depend on N. 

So, in the lowest order approximation (V = do + d 2 k 2 ) the spectrum of the short 
range interaction consists of one single mass M.. According to H, the low temperature 
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mass spectrum in k 4 approximation consists of two masses Ai + << M._ where Ai+ 
exactly coincides with M. presented above. 

Consider this point more carefully. Taking the case of v = 1 (N = 0) we present 
the zero temperature expression for T> up to the k 6 terms 

3/ 8\ 11/ 81\ 
tt 2 V(x) = 1 + 2ax ail 7 ) x 2 + — ail 7 I x 3 , (40) 



x 



2 V 3 7 18 V 11 



\ 2 2 2 \ 2 rfl 

a = 7in e X L , 7 = 7r n e X L £ . 



2-nnJi 2 ' 

The analysis of Q is based on the following k 4 approximation 

3 



7t 2 A(.t) = 1 + 2ax — - a ^— - 7^ x 2 



and the solution to the equation A(x) = leads to the conclusion about the existence 
of the second mass. 

In our case the k 4 approximation leads to the following factorization 

X?~ (k 2 +7W 2 )(k 2 -a 2 ), 

where the second factor is of a tachionic type. Situation can be clarified if the higher 
order terms are taken into account. For the typical values of the parameters we have 
7 ~ 0.0416 and one can easily verify that the equation 

(>n ~ 2 V'{x) = 11 (l - —7) x 2 - 18 ( 1 - - 7") x + 12 = 



a y ' V 11 7 V 3 

has no real solution, i.e. D[x) exhibits no extrema. This fact together with X>'(0) > 
implyies that T>'(x) > for all x. Moreover, we have T>(0) = 1 which means that 
defining equation T>(x) = has a single real solution x = Xq < 0. It is not difficult to 
guess that xq tends to —M 2 when we go to the lower approximation tt 2 T> = l+2ax. In 
other words, k 6 order corrections do not affect the lower order result: the short range 
interaction spectrum derived from (40) consists of a single mass defined by (39). The 
same time it is obvious, that the presented assertion depends on the approximation 
used and can be changed in higher orders. 

Constant parameters, corresponding to the long range components are given by 

s cs = {ng(o)ng(o) - nr(o)nS?(o) - £ ng(o)}, 



g u = ^n-(o), 

«0 



15 



and using (A. 6) - (A. 8) one can show that they vanish at T = for any u Q = N + 1. 

So, at finite temperatures the Coulomb forces carry the finite range character, 
while the charge -current and the current - current interactions include the short as 
well as the long range parts. When the temperature is decreased, integer number of 
the Landau levels are kept to be completely filled, and we see that the strength of the 
long range interaction decreases until vanishes at T = 0°!^, where the Maxwell field 
completely acquires the finite range character. 



Appendix 



Structure functions for a one type fermion system can be presented in the fol- 
lowing analytical form 



H E (x) = - A 1{x) _ _L e -*Ss(z,z) + I e"* A x JL So (x,y) 



, (A.l) 



y=x 



n cs(*) = ^ e 



d d 

dx~ Sl ( x ' v > ~ dx s ^ x ' v > + S ^ X ' V > 



y=x 



(A.2) 



n M (a?) = - tt" e 



d 2 



2n dxdy 
where I(x) and S a (x,y) are given by 



So(x,y) 



y=x 



(A.3) 



n=0a=l \ U °V • 



(A.4) 



S a (x,y) = @a(n)L n (x)L n (y)+ 

n=0 



Sa?i(« + «) ! 



e a (n,a)(x a + y Q )^(x)^(y). 



(A.5) 



Here are the generalized Laguerre polynomials, while 9 a are the temperature 

dependent constants 



aQ a (n, a) = 6 a {n) - 9 a (n + a), 



9 (n) = lim9 a (n,a), 

a— >U 



o (n) 



e^n) = a — + ep w 

OJn 



2 (n) = a £Pn+i, 
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9M 



1 / 

- UJ r , 



— )Pn+ o 



10. 



OJr. 



n+1 



+ 



p n+ i + cx£:(p n - p n+ i] 



Pn = Pn 0) + Pn (p), Pn = P„(p) - P n (p)- 

Structure functions which include the contributions from both sorts of fermions 
are denoted as 11^ (x), Hqs(x) and n^ t (x). For our purposes we consider the case of 
T = and at = 0. If the applied magnetic field B ext does not exceed the critical one, 
then at T = it is completely expelled from the sample (B h = 0), and according to 
(24) and (25) we have v\ = z/ 2 = v = u = N + 1 where N is any nonnegative integer. 
Consequently, only the nonvanishing terms in (20) and (21) in the zero temperature 
limit {(3 — oo) are 

Po(Pi) = ••• = Pat_i(pi) = P^(pi) = 1, 



P^(P2) = ••• = pjfM = pti+iM = 1, 

while all others, including the antiparticle contributions vanish exactly. These values 
of the Fermi distribution functions lead to the following zero temperature expressions 



ng(x) 



1 N ( X 3 1 

— V <^ -Ax + 3x 2 {2n + 1) {30n 2 + 30n + 11) }, (A.6) 

2?r n=0 I 9 J 

1 N ( St r 2 1 

= - E 1 " -7T ( 2n + !) + 7^ ( 30 ^ 2 + 30 ™ + n ) + 

W n=0 I 2 12 J 



3 iV 

+ — £(70n 3 + 105n 2 + 85n + 25), (A.7) 
727r n=0 



1 N ( 3r 1 

n^(x)= 2n+l- ^(2n 2 + 2n + l) + 

T n=0 1 2 j 



+ — E( 30rl3 + 30 ™ 2 + 32n + H) + 
127r n=0 



2 TV 



X 3 N 

+ E ( 35 ™ 4 + 70n3 + 120n2 + 85n + 25 )> ( A - 8 ) 

7 2?r n=0 

where the structure functions are expanded up to x 3 i.e. up to k 6 terms. Note, that 
the closed expressions for the nonrelativistic structure functions are presented in ||. 
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Figure 1. Free energies (conventional units) versus v for a = and T = 50°fT. 
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Figure 2. M [Gauss] versus B h [Gauss] for a = 0. 
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Figure 3. B h [Gauss] versus B cxt [Gauss] for a = 0. 
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Figure 4. Cricital magnetic field [Gauss] versus temperature [°K] for a — 0. 
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Figure 5. M [Gauss] versus B h [Gauss] for {-Kn e %/e) ■ a = 30 Gauss and T = 0°K. 
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Figure 6. M [Gauss] versus B h [Gauss] for {im e fi/e) ■ a = 30 Gauss and T = 40°^. 
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